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Abstract. We investigate spin effects in transport across fully interacting, finite size 
graphene armchair nanoribbons (ACNs) contacted to collinearly spin-polarized leads. 
In such systems, the presence of short-ranged Coulomb interaction between bulk states 
and states localized at the ribbon ends leads to novel spin-dependent phenomena. 
Specifically, the total spin of the low energy many-body states is conserved during 
tunneling but that of the bulk and end states is not. As a consequence, in the single- 
electron regime, dominated by Coulomb blockade phenomena, we find pronounced 
negative differential conductance features for ACNs contacted to parallel polarized 
leads. These features are however absent for an anti-parallel contact configuration, 
which in turn leads within a certain gate and bias voltage region to a negative tunneling 
magneto-resistance. Moreover, we analyze the changes in the transport characteristics 
under the influence of an external magnetic field. 

PACS numbers: 72.25.-b Spin-polarized transport; 72.80.Vp 
Electronic transport in graphene; 73.23.Hk Coulomb blockade; single- 
electron tunneling; 
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Figure 1. A graphene armchair nanostripc contacted by ferromagnetic leads. At 
the long sides, the lattice is terminated in armchair, at the small ends in zig-zag 
configuration. The length of a bond between two carbon atoms is do ~ 0.14 nm. We 
choose the orientation of the coordinate system such that the cc-axis points along the 
zig-zag ends, the y-axis along the long armchair edges of the stripe. 

1. Introduction 

Since the first successful separation of a one atom thick sheet of graphite by Novoselov 
and Geim and coworkers PQ, graphene, i.e. an isolated single sheet of graphite, has 
been a material attracting ever raising interest. Not only a great potential for appli- 
cations [21 [3], but also fundamental physics issues [I] arise from the linear dispersion 
relation in the electronic band structure of graphene, predicted about sixty years ago [5], 
as confirmed by various recent experimental findings [U El [TJ [8] . 

Increasing effort is presently put in the understanding of the electronic properties of 
graphene nanodevices, which can be obtained by etching or lithographic techniques and 
may achieve lateral dimensions of a few tenth of nm [91 [101 EI]- Studies on the effects 
of electron-electron interactions and confinement in transport across graphene nanode- 
vices have been carried out on single-dots [TQl EE2] and, only recently, on double-dot 
structures [131 E3] • A desirable goal is the fabrication of clearly defined geometries, and 
of particular interest for applications [2] could be narrow stripes of graphene, so-called 
carbon nanoribbons. Conductance quantization has been observed in 30nm wide rib- 
bons [TT], while an energy gap near the charge neutrality point scaling with the inverse 
ribbon width was reported in [9] and theoretically [151 HSj attributed to Coulomb inter- 
action effects. 

Crucial for the properties of a nanoribbon is the form of its ends. The two most 
regular possibilities are an armchair and a zig-zag end (see Fig. [T]). At present, the shape 
of the ends cannot be controlled, but there is ongoing progress in developing methods 
to fabricate stripes with clear edges, by scanning tunneling microscope lithography |17j . 
chemical synthesis |1_8] or by unrolling carbon nanotubes [19]. Moreover, there exist 
theoretical studies [201 [2TJ claiming that any narrow stripe should show either the be- 
havior of a zig-zag or of an armchair ribbon (ACN), where the names specify the form 
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of the long side edges. 

A peculiar property of the zig-zag edge is the existence of localized states [22] which were 
indeed observed experimentally by means of scanning tunneling microscopy j23j EU |25] . 
Due to their high degeneracy, flat-band ferromagnetism is expected from the Hubbard 
model, leading to a spin-polarized many particle ground state. It was suggested [261 EZ] 
to exploit this property for spintronic applications, where transport is governed by carri- 
ers in the oppositely polarized channels along the two long zig-zag edges. Recently both 
Hubbard and long-ranged Coulomb interaction effects have been analyzed [28] under the 
assumption of a filled valence and an empty conduction band (half-filling). There was 
a prediction of strong spin features in case of a low population of the localized mid-gap 
states. 

In contrast, in armchair ribbons the localized states are at the far apart zig-zag ends. As 
we showed recently [29], for narrow armchair ribbons short-ranged Coulomb interaction, 
i.e. local scattering events between either the two sublattices of graphene or between 
the extended bulk and the localized end states, gain increasing weight. They lead to an 
exchange coupling, inducing entanglement between the bulk and the end states. This 
has a decisive impact on eigenstates and transport properties. 

In this work we extend the investigations of Ref. [29] to the case in which the ACN 
is contacted to ferromagnetic leads and/or additionally subjected to a magnetic field. 
Various remarkable features due to the presence of the entangled end-bulk states, e.g. 
a negative tunneling magneto-resistance in a fully symmetric set-up, are discussed. 
The paper is structured as follows: In Sec. [2] and Sec. [31 we present our low energy 
theory for narrow ACNs. Readers not interested in technical details can, on the basis of 
the results for the ACN Hamiltonian summarized in Sec. I3.3[ directly start from Sec. HI 
where we shortly analyze a minimal set of states relevant for the explanation of features 
in transport across ACNs. Specifically, we show that scattering between end and bulk 
electrons causes an entanglement of states with the same total spin-projection S z , but 
different configurations of end and bulk spins. As a consequence, neither for the spins 
trapped at the ribbon ends, nor for the bulk electrons, S z is a good quantum number any 
more, which could be a handicap for proposed quantum information applications [2]. 
Although the formation of states symmetric or antisymmetric under the exchange of 
end spins remains without a special signature in the spectrum, it in fact leaves strong 
fingerprints in transport, as discussed in Ref. [29] for an unpolarized set-up as well 
as in Sec. |5]for collinear contact magnetizations. The most prominent feature in spin- 
dependent transport for collinear lead magnetizations is a negative tunneling magneto- 
resistance (TMR) within a narrow region along the edges of the Coulomb diamonds for 
even fillings, Sec. 15.21 Finally we explain how the transport characteristic is expected 
to change under application of an external magnetic field, both for non-magnetic and 
collinearly polarized contacts, Sec. [6j 
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Figure 2. The sublattice structure of the graphene honeycomb lattice. The underlying 
Bravais lattice is triangular, with a basis of two sites labelled p = ±. 



2. Noninteracting electrons in a finite size graphene armchair nanoribbon 

Being crucial for the subsequent analysis, we review here the end and bulk states of 
narrow, noninteracting ACNs. 

The carbon atoms in the graphene lattice are arranged in a hexagonal honeycomb lattice. 
There are two atoms per unit cell such that we can define two different sublattices p — ± 
as shown in Fig. [2J Hybridization of the 2s-orbital with the 2p x - and 2p y -orbitals leads 
to strong a-bonds in the lattice plane. The electrons in the remaining 2p z orbitals form 
7r-bands which determine the electronic properties at low energies. Characteristic for the 
structure of the 7r-band are the valence and conduction bands that touch at the corner 
points of the first Brillouin zone, also called Dirac points. Since there is one p 2 -electron 
per carbon atom, in isolated graphene the valence band is completely filled whereas 
the conduction band is empty. In the vicinity of the Dirac points the band structure 
exhibits a linear dispersion relation, Fig. |3K, resembling, up to a reduced propagation 
velocity of v = 8.1 ■ 10 5 m/s, the one of massless relativistic particles. From now on we 
focus on the region of linear dispersion. A description in terms of a Dirac equation for 
the p 2 -electrons capturing the essential features of the 7r-band close to the Dirac points 
can be obtained by a nearest neighbor tight binding calculation [30] . 
The six corner points of the first Brillouin zone can be decomposed into two subsets of 
equivalent Dirac points. As particular representatives we choose 



where ao ~ 0.14 nm is the nearest neighbor distance. Restricting the discussion to the 
vicinity of K±, the 7r-electrons are described by Bloch waves 



An 



k x ,F = ± 



(1) 



K F = F 






where f, R e M. 2 , Nl is the number of sites of the considered lattice, a = ± denotes the 
conduction and valence band, respectively, and Xr p (^) is the p z orbital on sublattice p 
at lattice site R. Finally, k is the wave vector relative to the Dirac point Kp. Defining 
the spinors 





Spin- dependent transport through interacting graphene armchair nanoribbons 



5 



it is found that they fulfill the Dirac equation 

hv (-Fa x K x + o-yKy) r] Fa (n x , k v ) = ae(n x , n y ) T} Fa (K x , (4) 

where 



e(k x , Ky) = kvJ k\ + K 2 (5) 



reflects the linear dispersion relation and o~ x , o~ y are the Pauli matrices. From (jlj) it 
follows that for K y ^ ±in x it holds the relation 

7] Fa + (K x , K y ) = -Ct— f 1 ' 2 = r)Fa-(K x , K y ). (6) 
V K x + K y 

A specific solution of Eq. ([6]) we will use in the remaining of the article is given by 

T)Fa-{Kx,Ky) = 1> VFa + {k x , Ky) = ~a- 



\f K x + K 



2.1. Boundary effects in narrow ribbons 



So far no boundary effects have been included. However, we wish to discuss the 
electronic properties of finite size carbon nanoribbons, which implies a wave function 
which vanishes all along the ribbon edges (open boundary conditions). The geometry we 
want to study is depicted in Fig. [TJ In particular we assume the long edges of the ribbon, 
along the y direction, in an armchair configuration, while the narrow terminations of 
the ribbon, along the x axis, have zig-zag character. We are interested in quasi one- 
dimensional ribbons and thus we restrict our discussion to stripes with a large aspect 
ratio L y ^> L x where L x and L y are the extensions in x and y direction, respectively. 

From Fig. |2]we can easily see that the two zig-zag ends each consist of atoms either 
sitting on sublattice p = + or p = — . We will use the convention that the "left" end 
at y = is formed by atoms living on sublattice p = — whereas on the other end (at 
y = L y ) we have atoms from sublattice p = + only. To fulfill the appropriate boundary 
condition for the zig-zag ends, we have to construct new wave functions ^(f, k) from 
linear combinations of ipFai?, k), in such a way that they vanish on the "missing" atoms 
at the ends, namely on sublattice p = + on the left end and p = — on the other end. A 
lengthy but straightforward calculation leads to 

(f> Fa (f, (K x , Ky)) 

= C zz (Fk x , K y ) [(p Fa (f, (k x , Ky)) - ip Fa (f, (k x , -Ky))} , (7a) 
with a normalization constant^ C zz (Fk x , K y ) £ C and the quantization condition [31] 



i2K v L y = FK x + IKy ^ 

%t\jy 

Additionally, the wave function we are looking for must also vanish at the armchair 
edges. In contrast to the zig-zag end, the terminating atoms where the wave function 

| The normalization constant is non-trivial as the functions given in Eq. (|2|) are non-orthogonal. To 
verify that there is the dependence on Fk x and k v , but not on a, Eq. ([2]) must be consulted in 
combination with Eqs. (TTJ) and (|6|). 
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is required to go to zero are from both sublattices. In order to build up suited linear 
combinations, we have to mix states of Eq. f JTaj) which belong to different Dirac points. 
Then, the resulting wave function (p a (r, k) vanishes on the lattice sites with R x = and 
R x = L x for 

(p a (r, (k x , Hy)) 

= C ac (n x , K y ) [{p K+a (r, (k x , K y )) - (p K _ a (f, (-k x , Ky))] , (8a) 
with C ac {n x , K y ) G C another normalization constant, and the quantization [3T] 

71 

K± ± k x = —n x , n x G Z. (86) 
L> x 

Due to the relation K± = ±K = ± 3% /| , the two conditions in Eq. (|86|) are equivalent. 
2.2. The eigenstates of metallic ACNs 

In total we obtain with Eqs. ( 17 ap and (18 qp the following expression for the eigenstates 
of noninteracting electrons in finite size ACNs, 

F,r=± 

where the generally complex number C(k x , k v ) = C zz (n Xl Ky)C ac (Kx, Ky) guarantees that 
V? a (r, (k x , K y )) is normalized to 1. 

We want to investigate now the solutions of the Dirac equation fulfilling our quan- 
tization conditions Eqs. ( l7fe|) and ()85|) . It has been shown [221 IMj that the presence of 
zig-zag ends leads to the formation of localized states characterized by a purely imagi- 
nary K y giving rise to an enhanced density of states around the Dirac energy. Though 
those states are localized at the zig-zag ends with an exponential decay in y direction, 
they are of decisive relevance for the transport properties of ACNs. 

At first, however, let us focus on the bulk states, where both k x and K y are real 
numbers. Since L y 3> L x , the quantization condition Eq. (j86j) leads to the formation 
of sub-bands assigned to different k x . For a metallic ACN, the sub-band must cross the 
Dirac points F = K±, which demands k x = 0, or equivalently n x = —K , n x G N, 
which follows from Eq. ( )8~6|) . The width L x depends on the number M of hexagons in a 
row parallel to x like L x = \/3ao(M + 1/2) (compare Fig. [TJ). Hence n x = |(4M + 2) 
(cf. Eq. (JTJ)) . Obviously, this means that the geometrical condition for having metallic 
ACNs with gap-less sub-bands reads Mmod3 = 1. Focussing on such ribbons, k cx k y 
as k x = 0, which corresponds to a cut through the Dirac cone as shown in Fig. [3|l. The 
corresponding states are characterized by the band index a and n y , see Eq. (l8~qj) . With 
n x = 0, Eq. (}76|) yields as allowed values of n y : 

Ky = y (n + , n G Z. (10) 




Figure 3. The dispersion relation of graphene, Eq. ([5]), for (a) real and (b) imaginary 
momenta K y . In both cases the boundary condition at the armchair edge mixes wave 
functions belonging to both Dirac points, (a) For solutions with real n y , in the low 
energy regime only sub-bands with k x = play a role due to the condition L x <C L y . 
The corresponding dispersion relation emerges thus from the intersection of the Dirac 
cone with the plane k x = 0. (b) For solutions with imaginary n y the cone opens along 
ik y with a slope determined by hv. As for the eigenstates it must hold n y = ±in x , the 
corresponding dispersion relation is obtained from the intersection of the cone with 
the two planes k x = ±ik y . There is only tangency along two straight lines within the 
k x -iky-p\ane, resulting in a dispersion which is identically zero. 



Since by definition, Eq. <fi a {r, (0, Ky)) = —(p a if, (0, —Ky)), we can further restrict 
for each a our analysis to either K y > or k v < 0. Thus it is an allowed choice to just 
consider states with sgn («;,,) = sgn(a), which we define as 

PKyir) ■= <fa= S gn( Kv ) (r, (0, Ky)) . 

Doing so, we select the positive slope of the two branches of the dispersion relation 
Fig |3K Bearing in mind the form of the graphene Bloch waves, Eq. (j2J), we can 
of course express the states ip b K in terms of the sublattice wave functions ^>f v := 

m L )- 1,2 T,R^ +K )\^{f), V 

Fpr 

where up to a complex prefactor the coefficients fp w are given by 

f F+r = rF, f F - r = iF. (11) 

Note that the index r here denotes right (r = +) and left (r = — ) moving waves [com- 
pare also to Fig. [3h] • 

Now we turn to the end states, emerging for purely imaginary K y [30], which are 
allowed by both the Dirac equation Eq. (j4j) and the quantization condition Eq. ( j76j) . 
In more detail, there exist two imaginary solutions for each k x > 1/L y , which holds 
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in ACNs for all k x = nw/L x , n G N. Besides, the relation L x <C L y causes that to a 
very good approximation K y = ±in x satisfies Eq. (j3J) and Eq. (l76j) . The corresponding 
dispersion relation is given in Fig. [3b. Notice that Eq. (EJ) is not applicable to n y = ±iK x 
(as explicitly exempted before), and instead the spinors fulfilling the Dirac equation Eq. 
(BJ are given by r)F ap (K x , ±i« z ) = S Pi ±f- Using this in Eq. (j2J) and following the steps 
leading to Eq. (Q, we obtain instead after straightforward insertions 

<Pa{r, {k x , ±ik x )) = ±C(k x , ±iK x )y*]Fp<p Fp (r, (Fn x ,ipn x )) . 

Fp 

The corresponding ACN eigenstates can be chosen such that they live on one sublattice 
p — ± only: 

F 

where C(k x ) is a normalization constant. It is evident that the decay length of <Pp Kx {r) 
from one of the ends to the interior of a specific ACN is L x / (n x 7r), which is always much 
shorter than the length of the ribbon. That is, one finds localized end states. 

From the dispersion relation Eq. (jHJ) it is easy to see that the energy of the end 
states is zero. Consequently, they will be unpopulated below half filling, but as soon as 
the Dirac point is reached, one electron will get trapped at each end (in an interacting 
system, Coulomb repulsion will hinder a second electron to enter). So we can conclude 
that at energies around the Dirac energy, not only the extended states with k x = 0, but 
also the localized end states can be of importance. 

2.3. Electron and Hamilton operator of the metallic ACN 

All in all, the appropriate operator describing an electron with spin a at position r reads 
in the low energy regime 

^(0= E E ( 12 ) 

Ky = (Z+0.5)n/L y p k x =Ntt/L x 

where c a k, d ap k are the annihilation operators for electrons of momentum k and spin a 
in the bulk or end states, respectively. The one- dimensional (ID) character of ACNs at 
low energies becomes evident by defining the slowly varying electron operators 

am--= E <,(o^ = ^= E ^ Ky ^y ( 13 ) 

such that we obtain 

VUr) = \A^/ 2 E fFpr<PFp(r)li>ra(y), (14) 
Fpr 

with ip Fp {r) ■= ¥Fp{r, K = (0, 0)). 



Spin- dependent transport through interacting graphene armchair nanoribbons 



9 



From the dispersion relations Fig. [3J it is easy to give the Hamilton operator of the 
noninteracting metallic ACN, 

H = hv^2K y cl Kv c aKy , (15) 

OKy 

with the Fermi velocity v = 8.1 x 10 5 m/s corresponding to the absolute value of the 
slopes of the linear branches in Fig. |3K There is no contribution of the end states as 
those have zero energy, see Fig. [3b. With the allowed values for K y , Eq. f JTOj) . there 
results obviously a level spacing 

e := hnv/Ly. (16) 
3. The interaction Hamiltonian 

This section is dedicated to the determination of the eigenstates of an interacting ACN. 
The cornerstones of the derivation and an analysis of the resulting spectrum were al- 
ready presented in our recent short manuscript [29] . We provide here merely a compact 
outline of the technical steps and refer the interested readers for details to the appen- 
dices. A discussion of the main features needed to understand the transport results is 
provided in Sec. [U 

In the following we concentrate both on interaction effects regarding the extended 
bulk states in ACNs as well as on correlations between end and bulk states. Ignoring 
exchange effects, the many-body end states can be spin- or edge- degenerate. Above half 
filling, within a reasonable energy range, both end states can be assumed to be populated 
with one single electron only. This is because their charging energy exceeds, due to the 
strong localization in position space, the charging energy of the extended states by far: 
A simple estimation modelling the Coulomb repulsion between two electrons localized 
at the same ribbon end via the Ohno potential given below, Eq. (Tl9|) . yields energies of 
the order of 0.1 eV for ribbon width around 10 nm. In contrast, typical charging energies 
for the bulk states of such ACNs range around 1 — 10 meV. The only relevant scattering 
processes between bulk and end states mediated by the Coulomb interaction are thus 
of the form 

^ = EE// d3rdV U{r-r')fc, § {r')Mr) 

crcr' p 

+ ^t(^%(r') U{f-r')^Ar')kp{ri) -(17) 

All other processes should be strongly suppressed for energy reasons. For the bulk-bulk 
interaction, our scattering potential is described by the expression 

V h . h = J J d 3 rd 3 r'$l(^l(r')U(r - F')^(0^(r ')■ (18) 

(7 (J 1 
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In both Eqs. (II 7p and (lisp , the function U(r — f') models the screened 3D Coulomb 
interaction potential. For our calculations we use the Ohno potential [32] , 

U(r- f) = U (l+( Uoelf 7j'\ ) ^ " , (19) 
V V V14.397[AeV]y J y J 

with Uq = 15 eV [33] and e ~ 1.4 — 2.4 [M] the dielectric constant of graphene. 
We proceed now with a discussion of the two different expressions. We start with the 
bulk-bulk processes, where the analysis follows largely the lines of an earlier work on 
interaction effects in metallic single wall carbon nanotubes (SWCNTs) [35] . 



3.1. Bulk-bulk interaction 

With the help of the reformulation of the 3D electron operator in terms of the 
ID operators ip ra {y), Eq. ( JT4|) . we obtain after integrating over the coordinates 
perpendicular to the ribbon axis an effectively ID expression for the interaction, 

Vb-b = ^2Yl l f ^^'^naC^^a'^O^^H^^O^a'^O^^), (20) 

{Ml ^' J J 

where ^{rm denotes the sum over all possible quadruples [h, I2, h, h], in the former 
case for the band index I = r. The spin-independent ID interaction potential U[ r ](y, y') 
reads 

L 2 

U[r](y,y') = -J- ^2 ^2 fhprifhp'r 2 fF3P'r 3 fF4 P r 4 
VP' {[F]} 

x ^dVdV^ Flp (rl^(fO^(f-fO^F 3P '(^)^F 4P (r), (21) 

with f indicating that the integration has to be performed over the coordinates 
perpendicular to y,y' (i.e. x,x',z,z'). Exploiting the explicit form, Eq. (ITT]) , of 
the coefficients fp pr we can easily identify those scattering processes which are indeed 
mediated by U[ r ](y,y'). In detail, 

U [r] (y, y') = 1 + ri ; 2r3r4 U^(y, y>) + ™ + r * r * U^(y, y'), (22) 

where we have defined the potentials 

L 2 

jjirara/inter^ ^ . = _± ^ ^ F 1 F 2 F 3 F 4 

P {[F]} 




x / / d rd r'^ p (^^± p (r / )C/'(r-r / )^Fs±p(r / )¥'J ;, tf (r)> 



describing interactions between electrons residing on the same/different sub lattices. 
From Eq. ( 122]) it is clear that a non- vanishing interaction potential can only be assigned 
to processes characterized by r^r^ = r 2 r 3 , and those are the forward (/)-, back (6)-, 
and umklapp (u)- scattering. Denoting the scattering type by Sj, the corresponding 
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quadruples are [I]g I= /± = [I,±I,±I,I], [I] b = [1,-1,1,-1] and [I] u = [1,1,-1,-1]. 
In total this means that we can rewrite Eq. fl20|) aq§| 

h-»= E E 

S r=u ,6,/± 5 CT =/± 

In the case of umklapp- and back-scattering with respect to r, the potential U[ r ](y, y') is 
proportional to the difference of the inter- and intra-lattice interaction potential. Since 
the latter potentials differ only on the length scale of the lattice spacing a$ ~ 0.14 nm, 
this means that in the case of S r = u,b the effective ID potential U\ r -\(y,y') can be 
considered as point-like. Introducing the coupling constant 



v 



l - J J dydy' (U intra (y,y')-U inter (y,y')) 



we can set in good approximation U[ r ](y,y') = u5(y — y') and write the short-ranged 
interaction processes as 

^/± b = \ E / d v4Uy)4U±*{y)$r±M1>-r<,{v), (23a) 

rcr 

K~- h = \ E/ dytUy)ti-M$-r-Mi>-r*(y)- (236) 

Since it is derived from a short-ranged interaction it scales inversely with the size of the 
underlying ribbon. We find typical values [36] of uL x /eq = 0.1 nm for a level spacing 
Eo, Eq. ( fl6l) . The process vanishes identically, because it involves the operator 

product i[>-ro-(y)'4 , -ra(y) = 0. Only the forward scattering processes Vf±*± are long 
ranged. 

Since easily diagonalizable by bosonization, it is convenient to identify the density- 
density processes among the relevant bulk-bulk interaction processes, such that Vb-b 
can be decomposed into Vb-b = V r p p~ b + V^ p ~ b , with the density-density and non-density- 
density parts given by 



-b-b _ T/b-b _| T/b-b j T/b-b _ T/b-b , T/b-b 

PP ~ / ± / ± + bf+ aIia V npp - V bf~ + V uf~ ' 



respectively. 



3.2. End-bulk interaction 

In ACNs, we additionally have to consider scattering between the electrons living in 
the bulk of the ribbon and the electrons trapped in the end state existing at both zig- 
zag terminations of the stripe. Below half-filling, the end states are unpopulated and 
thus all terms discussed in the following would be zero a priori. The range we want to 
concentrate on is the low energy regime above half-filling, where exactly one electron 
will permanently occupy each end state, so that we have in total two end electrons 
interacting with our bulk electrons. 

§ The spin quadruple [a, a 1 , a 1 , a] yields possible configurations [cr]f± = [a,±a,±a,a]. 
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We can certainly impose that the wave functions of the localized p z -orbitals have non- 
vanishing overlap for electrons on the same sublattice only. Moreover we demand that 
both end electron operators have to belong to the same end, and thus to the same 
sublattice, in order to give a nonzero contributior|Jj]. If we insert then into Eq. (fl7j) the 
decomposition Eq. ( fT4l) for the bulk electron operator and ipp a (r) = ( PpK, x i. r *)dc r p Kx 
for the end electron operator, and set in the coefficients from Eq. ffTTj) . we obtain: 

^ = t E E E E // d y <v 

Ki era' p rr' 

x (s rjr 4r'MU^(y,y%'p Kx - {rry^A f Ay')U:f{y,y%p Kx ) . (24) 

Thereby, symmetry arguments have lead to the demand r — r' for the part containing 
the interaction potential related to densities, 

U^ p (y,y') = J2 [ [ d 2 rd 2 r' FF'U{f-f') 

ppi J A- J _L 

x (^ + (r)^+(r) + ^_(r)^-(r)) \f e p Kx (r')\ 2 , 
which, again for symmetry reasons, fulfils further 

dy'U^ + (y,y') = [ dy'U^-(y,y') = t K *(y) « Uy) Mk x . 



We define a density-density part of the end-bulk interaction correspondingly as 

V e p ; h = L y Jdy t p {y) (^il{y)i ra {y)\ , (25) 

where we exploited that each end state is populated with exactly one electron. 
The second potential, 

U:f(y,y')= J2 [ [ d 2 rd 2 r'FF'U(f-f') 

ppi J A. J _L 

can be considered point-like due to the localization of the end states at yp=- = or 
j/p =+ = L y , and hence simplifies to 

Kf(y, y') = {5p,-5{y - 0)6(y' - 0) + 6 p>+ 6(y - L y )5{y' - L y )) 

with a short-range coupling constant which is independent of p for symmetry reasons, 

f n P := E // dVdVFF'f/(f-f0^(^^'+(O^(^)^(r)- 

ppi J J 

For an ACN of width L x ranging from 5 to 25 nm, one finds [3S] t^L ~ t np =: t, with 
tL x /eo ~ 0.55 nm, practically independent of k x . This means that the short-ranged end- 
bulk scattering is comparable in strength to the exchange interactions induced by the 

|| Operators acting on different ends would change the population of each end state, which is not 
allowed. 
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bipartite sublattice structure, and consequently we have to account for a non-negligible 
contribution 

Kp =-yE ^>+ rr ' + $Ul/p)&v(l/p) * E 4wU- ( 26 ) 

crcr'rr'p k x 

in the total end-bulk interaction potential V c ^ = V pp ~ b + V*~ p - 



5.5. Diagonalization of the ACN Hamiltonian 
express the total ACN Hamiltonian, 



We can diagonalize the Hamiltonian H +V pp h + V pp b by bosonization and subsequently 



#o := #o + V p h ~ h + V p c ~ h + V* pp h + V:; p \ (27) 



V v ' 



in the eigenbasis of Hq + V pp . A numerical diagonalization of the so constructed total 
Hamiltonian, however, yields reliable results only away from half-filling: as the eigenbasis 
needs to be truncated for the calculation, it is crucial that V npp is just a perturbation 
to Hq + Vpp i n the sense that it only mixes states close in energy to each other. In 
the direct vicinity of the Dirac points, the process V U JL breaks this demand, while it 
vanishes away from half-filling, as it will become clear in the course of this section. 

Diagonalization of the density- density part Diagonalization of Hq + V pp can be achieved 
by bosonization. The end result of the procedure on which more details are given in 
Appendix A is 



Ho + V PP = \ [W Af? + e tt c + [e - ^} ft J + E e jq a] q a jq . (28) 

The first three purely fermionic contributions in Eq. (128]) account for charging and shell 
filling. The last term counts the bosonic excitations of the system, created/annihilated 
by the operators dj g / hj q . The two channels j = c, s are associated to charge (c) and 
spin (s) excitations. The excitation energies are 



£ jq — \ X jq A ]q ' 

u 

with X jq = n q [S jtC W q + e ] , A jq = n q 5 jtC W q - - . (29) 

The energies of the charge channel are dominated by the long-ranged interactions via 
the coefficients 

W * = T2 f LV dy' (U mtra (y,y') + U mter (y,y')) C os(qy)cos(qy'). (30) 

L y JO JO 

Finally, we can give the eigenbasis of Hq + V pp ~ h in terms of states 

N,a e ,m)= lHl( a l) mj9 Mv|^<? e ,0) , (31) 

j=c,s q>0 
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where \N,a e ,0) has no bosonic excitation. The fermionic configuration N = (iV-j-, 7VJ 
defines the number of electrons in each spin band. The occupation of the end states 
is determined by a e = (o~+,crL), where '— ' relates to y_ = and '+' to y + = L y as 
before. Below half filling the end states are empty, such that there is only one possible 
configuration: — = o~ e _. Above half filling, exactly one electron occupies each 
end state and thus a\,a e _ G {til}- Finally, m = (m s ,rh c ), with rrijq = {fhj) q contain- 
ing the information how many bosonic excitations are present in level q for mode j = c, s. 



Non- density- density interaction In the following we use the states from Eq. (l3Tj) as 
basis to examine the effect of V npp . For this purpose we evaluate the matrix elements 
of the potentials V^~ p h and V° pp , using the bosonization identity for the ID electron 

does not conserve the quantity m, while it must neither mix states 



operators. 

Generally, V^~ h 

with different electron configurations N, nor with different end spin configurations 

a e 

rb-b „„j T^b— b „„„j_„; J ;„ T>b— b 



the Coulomb interaction between bulk electrons cannot change the quantity 
= |(A^- — N±), and it cannot touch the end states. Further, we already know 
that both processes V^-y and V; b " b contained in V^~ p h are effectively local interactions, 



see Eqs. (123 a\i and (1236ft . such that the matrix elements of the non-diagonal bulk-bulk 
interaction scale with the exchange-coupling parameter u, 



Na e m K b " b 



npp 



N'a e 'm 



u 



2£ w N,N' w ™c,m' c <- 



X 



e 



4 sin (7ry/L y ) 



, m sq , m 



sq i 



(32) 



The derivation of this expression is given in |Appendix B[ as well as the definitions 
of F(X,m,m') and A^jr ,(2/), Eq. (1B.14|) and Eq. (1B.12|) . respectively. In fact it turns 
out that in comparison to the end-bulk non-diagonal interaction, the bulk-bulk non- 
diagonal interaction has only minor impact on spectrum and transport properties of 
narrow ACNs. 

For the non-diagonal end-bulk interaction, we find, as explicitly evaluated in 



Appendix C 



(N a e = (o+,o-_) m 



V. 



e-b 
npp 



N' a e ' = (a' + ,o J _) m! 



.,m' c ^cr p , — cr' p dap , a' p 



N> 



sgn(c7- p ) J]F Usgn((/ p )J2/ 



'n q ,m sg ,m sq 



(33) 



1 

The action of this scattering process is to flip the localized spin at the p end, o~ p — > a p = 
—o~' while at the same time a bulk spin must be inverted to preserve the spin-projection 

S z : N' a , ->• N a > = N' a , + 1 , N'_ a , N^ a , = N'_ a , - 1. The localized spin at the other 
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end, i.e., the p = —p end, must be conserved: a p = a' p . As a result, the degeneracy 
between the lowest lying states of Hq + V pp is removed; moreover, the spin-charge sepa- 
ration gets smeared [29]. 

4. Minimal model for the lowest lying states 

The low energy spectrum resulting from the numerical diagonalization of H Q , Eq. ( |2"7|) . 
is discussed in Ref. [29]. Here to explain how the lowest-lying states in the truncated 
eigenbasis Eq. ( )3T]) transform under the influence of V npp it is sufficient to restrict to a 
minimal set of states: For the even fillings, N c = 2n, n £ N, we have to take into account 
twelve states, allowing up to one fermionic excitation. The reason is that without an 
unpaired bulk spin no mixing can take placed- For the odd fillings, N c = 2n + 1, it 
is enough to include the eightfold degenerate ground state of Hq + V pp . To preserve 
lucidity, bosonically excited states are left out from our analysis, as they do not change 
qualitatively the mechanisms behind the observed effects. The N and m c conserving 
bulk-bulk scattering V^~ p h can also be disregarded: besides for slight shifts in energy, 
forming linear combinations of states differing just in m s but identical in N, rh c has 
not much impact. The restriction to the minimal model applies only for the subsequent 
analytics. Concerning the numerical calculations shown in Sec. an energy cutoff of 
1.9eo above the ground states was used, including every energetically allowed bosonic 
or fermionic excitation. 

4-1. Even electron fillings 

For even electron fillings, N c = 2n, we want to restrict to the following subset of the 
states described by Eq. (J3lj) : 

N = (N t ,N i )e{(n,n),(n±l,nTl)} n £ N, 
<? e = «,ai) £ {(t,t),(t4),a,t),(U)}, 
rh = (rh c ,rh s ) = (0,0). 

To abbreviate our notation we introduce 

[ \ot,n,°%) N t = N x = n, 
iV,<T e ,0):=<^ k,tt,<> N t = n + 1, N i = n-1, 
[ \a e _, U,o+) N t = n-1, N^n + 1. 

Notice that the second and the third case describes fermionically excited states. Building 
now all possible combinations for our minimal set of states for even fillings, we get the 
following set of possibilities: 

For transport, this is fatal, as introducing the end spins degree of freedom a priori yields four 
identical, completely decoupled channels. This makes the kinetic equations ill-defined. Only the end- 
bulk interaction induced mixing guarantees a unique solution for the transport problem. 
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S z = 0: four states \a), \b), \c±) 

ft,n,l)=:\a), \i,UA)=-\b), 
|t,|4,t)=:|c+), \i,n,i)=-\c-), 
S z = ±h: six states \d aa ), \f aa ), \g aa ) (a G {|,4}) 

tt,UA)=--\dn), \l,U,l)=--\d u ), 
It,tt4)=:|/n>, lt,W4)=:|/u>, 

|l,tt,t)=:|0tt>, \i,UA)=-\9u), 
S z = ±2h: two states \a, a a, a) (a G {1\4}) 

|t,tt,t>, \l,U,D- 



(0) 



E®, while 

(/) 

2n • 



There are four degenerate ground states \a), \b), \d aa ) with energy E 
the remaining singly fermionic excited states \f aa )i Isw), \°~, crcr , °~) have an energy E. 
The end-bulk interaction can only mix states with same spin-projection S z . To the 
highest values, S z = ±2h, it belongs only one state and thus no mixing occurs. In 
contrast, the three states \d aa ), |/ CTCT ), |gw) with S z = sgn(cr)ft get coupled to each 
other. Also the four states \a), \b), \c±) with S z = can in principle transform into 
one another under the influence of end-bulk scattering. With help of Eq. ( 1331) we can 
set up the corresponding blocks of the full Hamiltonian H Q in the truncated eigenbasis 
(assuming n evero) 



#0 



N c =2n,S z =0 



\a) 


\b) 


\c + ) 


Y-) 






( E {0) 





-t 


-t 


\ 


\a) 







+t 


+t 




\b) 


-t 


+t 


E (f) 







|c+ 


{ -t 


+t 





E U) 

- C/ 2n 


) 


c_ 



N c =2n,S z =±H 



) \Ua) \9ao) 

( Eg +t -t \ 







E. 



if) 

2n 



\daa) 
\f«a) 



Diagonalization leads to the eigenstates and eigenenergies listed in Tab. [TJ We employ 
there the abbreviation 



6aa'(7) 




aEg + a' 



E if) - E 



(0) 



r 



(34a) 



with a, a' G {±1}- Obviously, £++(7) > £+-(7), £_ + (7) > £ (7), and as in our 



context 7 — t, and hence 7 <C E^J — E^l holds, we can rely on the relations 



(0) 



£++(7) n E%} , e+-(7)^^ u J, e-+(7)»t, e-(7)«t. 



,(0) 



(346) 



For odd n the sign of the off-diagonal entries must be inverted. 
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N c = 2n 

ENERGY : EIGENSTATE (NOT NORMALIZED) ABBR. SPIN S z [H] 

£ + -(4t)^0 : i 2v|_(|a)-|6))-(| C+ ) + |c_)) | gl > | 

£ + _(2V2*)~0 . ^y 2t) \ daa ) + {\f aa ) -\ gaa )) \g 2 ) aa # sgn(a) 

: y/l(\a) + \b)) |gs) ■ 

k+(4t)« eo : (|a> - |6» - (|c + ) + |c_)) | ei ) ■ 

Ein^eo ■ jsdM + lQaa)) • Sgn(<7) 

Eg?«e : y^(|c + ) + |c_)) ■ 

#2? ~ £ o : k, o-o", a") ♦ 2sgn((x) 

e ++ (2y^t)^ g0 : g _ + 2 ^ t) |^) + (|/ CTCT )-|^)) |e 2 ) gg # sgn(a) 



Table 1. The minimal model consists of twelve different states for even electron fillings 
of the ACN (a e {t)4})- We have set £^ = 0. Importantly, fermionically excited 
states mix with non excited ones. States relevant for explanations in the main text are 
marked with text labels. The color and shape of the symbols encodes the nature of the 
state, classified by spin-projection and behavior under exchange of the end spins. The 
boxes indicate S z = 0, related colors are red (symmetric component \a) + and blue 
(antisymmetric combination \a) — \b)). For the disks \S Z \ = H, colors olive (symmetric 
components) and green (antisymmetric: \f aa ) — \g a a))- Grey diamonds: |iS^| = 2h. 

The resulting energy landscape is sketched on the left side of Fig. HJ where we used 
differently colored and shaped symbols to indicate the composition of states. In our 
simple model, we find then from Tab. [1] that the interaction has hardly lifted the 
degeneracies between the various states. It can be verified with Eq. (134 aft that there is 
a slight splitting in the ground states, such that their energy grows from |gx) to |g 3 ). 
Also the excited states, of which only the two labelled ones turn out to be important for 
later explanations, are listed increasing in energy. Crucial is the mixing of states with 
different bulk and end configurations. We can single out linear combinations which are 
either symmetric or antisymmetric under exchange of the end spins. For the features 
we will observe in transport, the decisive entanglement is the one between the four 
S z = states \a), \b), |c+), |c_), leading to two states containing the antisymmetric 
combination \a) — \b): A ground state |gi) with small contribution of the fermionically 
excited states \c + ) + |c_), and an excited state |ei) where those dominate [as found with 

Eq. (Jhsd;. 
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2n 



Energy 



E W> + £ 



E 2n ^ 



ei 



IgO |g*> Iga) 



|e 2 > 



2n+l 

Energy 



State symbol colors: 
Go-Sow component 
c^/d aa /^ a +g aa component 
a+to component 

a-b component 
|Sz| maximal 

State symbol shapes: 
|Sz| maximal O 
|Sz| minimal □ 
in between O 



Figure 4. Energy landscape for a minimal set of lowest ACN eigenstates for even 
(N c = 2n) and odd (N c = 2n + 1) electron fillings, in accordance with Tabs. [T] [2] 
To visualize the relevant contributions in the composition of the eigenstates, different 
colored and shaped symbols were used. The states relevant for later considerations are 
labelled. 



4-2. Odd electron fillings 

Here, due to the fact that with N c — 2n+ 1 we necessarily always have an unpaired spin, 
it is sufficient to consider merely the ground states with energy E$ = E^ +l emerging 
from Eq. (1311): 

N = (N^N±) e{(n+l,n),(n,n+l)} neN, 
a e = K,a!„) G {(t,t),(t,l),(l,t),(l,l)}, 
rh = (rh c ,m s ) = (0,0). 

Again we abbreviate our notation and introduce 

AT B e o\ •=/ \ a -^> a +) N t = n+l,N± = n, 
' ' / • \ \<J±A,°%) N t = n, N i = n + 1. 

We get the following set of possibilities: 

• S z = ±h/2: six states \a a ), \b a ), \c a ) (a =f-H- a =|, a =|<H- a =|), 

\t,cr,D -.\a a ), ||, a, t) =: |&<r), k a, a) =: |c CT ), 

• S z = ±3/1/2; two states |ct, a, a) (a e {t,l}) 

|t,t,t), 11,1,1)- 
In the case of S z = ±3h/2 there is only one state each. 

For the three states with S z = ±h~/2, from Eqs. (1281) and (133]) the following mixing 
matrix is found (still n is assumed even): 
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N c = 2n + 1 



ENERGY 


EIGENSTATE (NORMALIZED) 


ABBR. 




SPIN S z [h] 


-V2t 


§(K> + |6«r» + \f\K) 


|tl) CT 


■ 


sgn((j)/2 









■ 


sgn(fi)/2 





1 <r, o-, a) 




♦ 


3sgn(<r)/2 


+V2t 


f(k> + \b*)) - \J\K) 


|t 3 ) CT 


■ 


sgn(a)/2 



Table 2. Lowest lying eigenstates of an ACN filled with an odd number of electrons 
(a G {^,^})- Due to spin-degeneracy, the total number of possible states is eight. The 
two states with \S Z \ ~ ih/2 are marked by grey diamonds. Red and olive boxes stand 
for the symmetric components \a a ) + \b a ) and \c a ), respectively. Blue boxes label the 
antisymmetric combinations \a a ) — \ b a ) . Notice that all states behave purely symmetric 
or antisymmetric with respect to end spin exchange. 



N c =2n+l,S z =±K/2 





\a a ) 


K) 


\Ca) 






/ 


E {0) 





-t 


\ 









p(o) 

rj 2n+l 


-t 




K) 


V 


-t 


-t 


E (0) 
nj 2n+l 


J 


\ c a) 



The matrix is easily diagonalizable and yields eigenstates according to Tab. [2] at three 
distinct eigenenergies (compare also to Fig. HJ right). Notice that for the odd filling all 
emerging states are purely symmetric or antisymmetric under exchange of the two end 
spins. Thereby, the symmetric states |ti) CT and ^3)^ essentially have the same tunneling 
properties, because they only differ by the sign in front of \c a ). It is of crucial importance 
that the state ^2)0- is formed by the antisymmetric combination \a a ) — \b a ). Comparing 
the definition of the 2n states \a), \b) and \a a ), \b a ), we see that from the 2n ground 
state |g 3 ), a tunneling event can exclusively lead to one of the 2n+l states |ti) CT or ^3)0-. 
Via their \c a ) components, these connect to \c + ) + |c_) as well as to \d acr ), and thus to 
all the other labelled 2n states from Tab. (TJ but the link to |gi), 162)0-0- is weak due to 
Eq. ( 13461) . This will be the key ingredient for the explanation of the stability diagrams 
in Sec. El 



5. Spin-dependent transport across quantum-dot ACNs 

Looking solely to the spectrum, there is no demand for distinguishing between symmetry 
or antisymmetry of a certain state under the exchange of end spins. In transport, how- 
ever, this property turns out to lead to tremendous effects. The case of an unpolarized 
set-up was discussed in Ref. [29], where it was found that end-bulk entanglement leads 
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to pronounced negative differential conductance (NDC) lines occurring for a completely 
symmetric setup. In this work the focus is on spin-dependent transport for collinear lead 
magnetizations. Strikingly, all the NDC features observed for the unpolarized set-up 
vanish for anti-parallel contact magnetization while they persist for the parallel case. 
As a consequence we predict negative tunneling magneto-resistance (TMR) within a 
narrow region along the edges of the Coulomb diamonds for even fillings. 

Unless specified differently, we have employed the following parameters for all 
viewed plots: 



Energy cutoff 


E 


1.9e , 


Thermal energy 


k B T 


0.01 meV, 


Charging energy^ 




2.31 meV, 


Ribbon length 


Ly 


572 nm, 


Level spacing 


eo 


2.93 meV, 


Ribbon width 


L x 


7.8 nm, 


Bulk-bulk exchange 


u 


0.036 meV, 


End-bulk exchange 


t 


0.21 meV. 


Polarization strength 


P 


0.8 



The values of charging energy, bulk-bulk and end-bulk exchange-coupling were numeri- 
cally verified for ribbon widths ranging from 5 — 20 nm. 

Throughout this work we assume that the coupling between the electronic 
reservoirs, i.e. the contacts, and the ACN is weak. Under such condition, the total 
system is described by 

H = H e + H j ea( j s + H T — eaV gatc Af c , 

with the ACN-Hamiltonian H Q given in Eq. (1271) . Further, the contacts are described 
by H\ eac [ s = Yliq ^2a( e q~ ^claqCio-q, with C\ aq annihilating an electron in lead / of kinetic 
energy e q . The chemical potential \x\ differs for the left and right contact by eVu&s, with 
Vbias the applied bias voltage. Next, Ht = J2ia ffir (Tiffifyffify^r) + h.c. j describes 
tunneling between ACN and contacts, where Ti(r) is the in general position dependent 
tunneling coupling and ip a (r) the ACN bulk electron operator as given in Eq. (fH|) . The 
lead electron operator is <t>i ai {r) = Tlq^iqif)^^ with <fi q (r) denoting the wave function 
of the contacts. Finally, the potential term describes the influence of a capacitively 
applied gate voltage (0 < a < 1). Due to the condition that the coupling between ACN 
and the contacts is weak, we can calculate the stationary current by solving a master 
equation for the reduced density matrix to second order in the tunneling coupling. As 
this is a standard procedure, we refer to previous works [33, EHl [39] for details about 
the method. 

* A dielectric constant e = 1.4 was assumed. 
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Besides the energy spectrum, the other system specific input required for transport 
calculations are the tunneling matrix elements of the ACN bulk electron operator, 

1 



iVa e m|^ CT (f) 

X ^ fFpr^Fpi 
Fpr 



N'a e 'm! 

r e L v 



2 



8-„ 



(-1)' 



\)ry 



q>0 j=s,c 



iq (y),m jqi m' jq ) 



with the function F(A,m, m') and the parameter X J r q T (y) as given in Appendix B in Eq. 
(IB. 141) respectively Eq. (IB. 111) . We omit here the calculation of this identity, because 
the fermionic contribution follows straightforwardly from Eqs. (IB.lj) - flB.3l) . while the 
bosonic contribution emerges in the same manner as for the more complicated matrix 
elements involving more than one electron operator evaluated in Appendix B Moreover, 
the detailed derivation of the corresponding expression for SWCNTs can be found in 

Ref. na. 



5.1. Collinearly spin-polarized transport without magnetic field 

In Fig.[5j left, (a) and (b) we show the stability diagrams obtained for an ACN coupled to 
leads polarized in parallel and in anti-parallel, respectively. Due to electron-hole sym- 
metry around the In filling, the transport characteristics are mirror-symmetric with 
respect to the central diamond. As anticipated at the beginning of this section, end 
states leave various signatures in the parallel case, Fig. El left, (a) which are absent in 
the anti-parallel configuration, Fig. [5], left, (b). In the following we give an explanation 
for all effects indicated in Fig. [5], left, (a) based on the minimal set of states discussed in 
Sec. 14.11 Tabs. HJ [2j Let us shortly recall their properties and alongside explain in which 
points we have to expect discrepancy with respect to the full set of real eigenstates 
which was employed for all calculations: 

Firstly, the eight states from Tab. El |ti) CT , \t 2 ) a , ^3)0 an d \cr, a, a), will in the following 
frequently be called the lowest lying In + 1 states. They occur at only three distinct 
energies ±y/2t, 0.0. Inclusion of excitations within an energy cutoff of 1.9 Eq slightly 
lifts the degeneracy of \t 2 ) a an d |cr, er, a), and introduces eight high lying excited states 
which are almost degenerate. Secondly, for what concerns the even fillings, we refer to 
|gi) , \g2) aa, \g3) as 2n ground states. The fact that they are almost, but not perfectly 
degenerate, and that their energy grows from |gi) to |g3) , is not changed upon the in- 
clusion of further excitations and plays some role in the following. Moreover, mixing 
between the states from Tab. [Hand bosonically excited 2n states takes place in general, 
but actually preserves the types of linear combinations occurring in Tab. [H which is 
the relevant point for our explanations. In summary, the main effect of inclusion of 
excitations within an energy cutoff of 1.9 Eq is the lifting of the degeneracies among the 
excited 2n states. In fact, the lowest lying excited state will be of the same nature as 
|ei), and the separation to the state corresponding to 162)0-0- exceeds y/2t, thus it is well 
resolvable. With this additional information to Tabs. [U [2] we can now start to explain 
the features marked by the different arrows in Fig. [5] (left). 
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(a) Parallel polarized leads, P=O.S 



2.5 5.0 
Gate voltage eaV^ E [meV] 



N c =2n+2 



Diff. cond. 

d//c%„ 
[10-' eVfi] 



0.4 
0.0 



Population of |gz) Occupation 
probability 

' 1.0 



0.6 0.8 1.0 1.2 
Gate voltage eaV g , K [meV] 



0.8 
0.6 
0.4 
0.2 
0.0 



(b) Anti-parallel polarized leads, P=0.8 




2.5 5.0 

Gate voltage ea^ ( [meV] 



Population of |g 2 > 0ccupation 

probability 
1.0 



0.6 0.8 
Gate voltage 



1.0 1.2 

mV [meV] 



0.8 
0.6 
0.4 
0.2 
0.0 



Figure 5. (a) Left: Differential conductance for an ACN transport setup with 
parallel lead polarizations. Strong negative differential conductance (NDC) arises if 
the external voltages are adjusted such that a transition from Iti)^ to |ci) is allowed, 
while a transition from |gi) to |t2)cr is forbidden, (a) Right: Occupation probability of 
the trapping state |gi) around the region exhibiting various NDC features. Notice that 
no numerically stable data can be obtained inside the Coulomb diamond, (b) Left: 
Differential conductance for an ACN quantum-dot connected to anti-parallel polarized 
leads. The number of visible transition lines is strongly decreased as compared to 
the parallel case Fig. [5ji (left). Further, all NDC features have vanished, (b) Right: 
Occupation probability of the state |gi) for the same bias and gate range as in Fig. 
[5ja, (right). For an anti-parallel contact configuration, the population of the state is 
strongly decreased. 



The dashed red arrow points towards a triple of three parallel lines which are 
split by \/2t. Those mark transitions from the 2n ground states to the In + 1 lowest 
lying states. Hereby, the antisymmetric state ^2)0-, associated to the second line of the 
triple, is special, because it is the only one strongly connected to the 2n state |gi). The 
first line of the triple is the |g 3 ) — >\ti) a ground state transition line. 

The blue dotted arrow marks the lines around the tip of the Coulomb dia- 
mond. Those appear in four clearly distinct positions, separated by about \/2t . The 
lower lying triple of lines arises from transitions of the lowest lying In + 1 states to the 
2n + 2 ground states. By coincidence of parameters, the highest line, which is split, 
follows also in a distance of about \f2t and marks transitions from 2n ground states to 
the aforementioned higher lying excited 2n + 1 states arising upon inclusion of bosonic 
excitations. 
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Figure 6. Schematic explaining the mechanisms causing the NDC features (I), (II) 
and (III) in Fig. [5] (left). Only states and transitions relevant for the NDCs are drawn. 
The crucial transition is marked by a big arrow head. (I) / (II) Opening of the channel 
|ti)<j— >|e2)(T<T, respectively |ti) CT — >|ei), leads to a decay into the trapping state |gi), 
depleting the transport channel |g3)<-^|ti)o-. (Ill) Opening of the channel |t2)cr — H e i) 
depletes the transport channel | gi )•(->■ |t2)o- . 



The solid green arrow highlights the negative differential conductance features 
(I), (II) and (HI). The former two originate from trapping in the state |gi), while the 
latter occurs due to depletion of the transport channel \t2) a ^\Si) ■ 
The mechanisms work as follows: 

The NDCs (I) and (II) mark the opening of the 2n + 1 — > 2n back-transition channels 
|ti)o— H e 2)o-<7 and |ti)o-— Hei), respectively. The situation is sketched in Fig. [6j Once 
they get populated, from both of these excited 2n states the system can decay into any 
of the lowest lying 2n + 1 states, and in particular there is a chance to populate the 
antisymmetric state ^2)0-- This state is strongly connected to the 2n ground state |gi), 
which contains a large contribution of the antisymmetric combination \a) — \b). But 
in the region where the NDCs occurs, the forward channel |gi)— >-|t2) CT is not yet within 
the bias window such that |gi) serves as a trapping state. Fig. [5^ (right) confirms this 
explanation: the population of the state |gi) is strongly enhanced in the concerned re- 
gion where the back-transitions |ti) CT — >\e2)ao- and |ti)o-— >• |ei) can take place, while the 
forward transition |gi)— )-|t 2 ) fJ is still forbidden. 

NDC (III) belongs to the back-transition |t2)o- — which is a weak channel because 
^2)0- is a purely antisymmetric state, while the antisymmetric contribution in |ei) is 
rather small. From time to time, nevertheless the transition will take place, and once 
it happens the system is unlikely to fall back to |t 2 )o-, but will rather change to a sym- 
metric 2n + 1 state. Thus the state |t 2 )o- is depleted, and with it the transport channel 
^2)0-^ |gi) , which leads to NDC. The statement can also be verified from the plot of the 
occupation probability for |gx) , Fig. [5] (right): a pronounced dark region of decreased 
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population follows upon the NDC transition. 

Major changes in the stability diagram of the ACN are observed for anti-parallel 
contact configuration, Fig. [5b (left): Compared to Fig.[5k (left), various transitions lines 
are suppressed and the NDC features have vanished. The reason is that an anti-parallel 
contact configuration as drawn in Fig. [1] favors ^2)44 as 2n ground state, because in- 
tunneling of 4- - electrons and subsequent out-tunneling of 7 - electrons is preferred. As 
a consequence, all transport channels related to |gi) and |g3) are of minor relevance, 
which weakens various transition lines and in particular destroys the NDCs effects: in 
the anti-parallel configuration, the occupation of the trapping state |gi) is significantly 
lowered, as seen in Fig. [5b (right). 

In detail, starting from the 2n ground state |g3)= ^(\ a ) + |^)) ; in-tunneling of a ma- 
jority (4 -) electron from the source takes the system to (ti)^ = -^(\ a i) + I&4)) + |ci). 
Via the |cj,) (= 1 1, f, 4)) - component of this state, it is possible to tunnel out with 
a majority (7 -) electron of the drain, yielding a transition to (§2)44.- Similarly, also 
starting from \gx) in-tunneling of a 4 - and subsequent out-tunneling of an 7 - electron 
changes the In ground state to ^2)44- Depending on the bias voltage, transport is either 
carried by 7 - electron via the ground state channel ^2)44 -H- (ti)^, or by 4 - electrons 
via (g2)44 -H- 1 4> 4> 4-), where |4, 1, 4) forms a blocking state unless a back-transition to 
the 2n excited state |4, 44, 4) is energetically allowed. 



5.2. Tunneling magneto-resistance 

In Fig. [7] we have plotted the tunneling magneto- resistance (TMR), 

TMR := Ipa ~ Iap j 
Iap 

which is a measure for the ratio of the current in the parallel configuration, Ipa to 
the current in the anti-parallel configuration, Ipa- Along the edge of the 2n Coulomb 
diamond, the TMR acquires a negative value, i.e. Iap exceeds Ipa- This is unusual: for 
lowest order calculations without Zeeman splitting between the spin species typically 
strictly positive TMR is observed jHJ [381 EH]- For the ACN, however, the effect 
originates from a reduced feeding of the |gi) trapping state. This statement can be 
confirmed by comparing its occupation probability for the parallel, Fig. [5^ (right), and 
anti-parallel polarized case, Fig. [5b (right), in the concerned region. 

Namely, |gi) lies slightly lower in energy than |g3) -for the values we chose, the 
energy difference amounts, according to Tab. [I] and Eq. (134 ap . to about 6kpT. Hence 
it can serve as a perfect trapping state within a narrow region along the edge of the 2n 
Coulomb diamond: here, the bias is high enough to allow the ground state transition 
|g3) — ►|ti) (7 , but not |gi)— >|ti)a-. Though the latter channel is weak in any case, it 
nevertheless provides a nonzero escape rate from |gi). That is why in Fig. [5^ (right), 
in the region where the NDC mechanism Fig. [6] (II) can populate |gi), the occupation 
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0.0 2.5 5.0 7.5 

Gate voltage eaV aM [meV] 

Figure 7. tunneling magneto-resistance Ipa/Iap — 1, where Ipa/ap denotes the 
current for parallcl/anti-parallcl polarized leads. We observe a negative value, i.e. 
Ipa < Iap, along the edges of the 2n/2n + 2 Coulomb diamonds as soon as the channel 
|ti)o-— Hei) has opened. The reason is the decreased population of the |gi) state in the 
anti-parallel case [Fig. [SJd (right)] as compared to the unpolarized or parallel case [Fig. 
Ek (right)]. 



probability approaches 1 only straight along the edge of the Coulomb diamond, and a 
value of 0.6 — 0.8 further apart from it. In contrast, we observe no comparable increase 
of the |gi) population in Fig. [5b (right), because for the anti-parallel configuration, as 
explained above, the transition channels involved in the NDC mechanisms are strongly 
disfavored compared to ^2)44 -H- (ti)^. For this reason, no trapping occurs and Iap can 
exceed I pa, leading to the negative TMR. 

6. Magnetic field sweep 

Finally we study, both for non-magnetic and collinearly polarized contacts (see Fig. 
[I]), the transport behavior under the influence of an external magnetic field. In its 
presence, formerly degenerate states with different spin projections S z components 
become Zeeman split. This means, at a fixed gate voltage one half of the transitions 
occur at a higher, one half at a lower bias compared to the situation without magnetic 

field. In detail, simple thoughts can confirm that the forward transitions involving j" - 

t 

electrons, i.e. 2n — > 2n+ 1, as well as all backward transitions mediated by J, - electrons, 
i.e. 2n + 1 — > 2n, are lowered in bias, while processes 2n — y 2n + 1 and 2n + 1 — > 2n are 
raised. 

For three distinct values of the Zeeman splitting, Fig. [8] shows stability diagram 
zooming on the region between the 2n and 2n + 1 Coulomb blockade diamonds. Those 
plots complement Fig. Efe., where we show the differential conductance versus bias and 
Zeeman splitting, at a fixed gate voltage eaV gSLte ~ 2.0 meV (marked in Fig. [8] with the 
dashed line). In turn, the three values of the Zeeman splitting considered in Fig. Elare 
marked in Fig. |9^i by dashed white lines. 
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Figure 8. Differential conductance in the region between the 2n and 2n + l Coulomb 
blockade diamonds for an ACN quantum-dot at three different values of an external 
magnetic field. The contacts are assumed non-magnetic, i.e. unpolarized. The dashed 
white lines mark the gate voltage at which the plots Fig. |9]are taken. 



At first we focus on Fig. [9^- In Fig. [9b the case of polarized leads is considered. 
A small Zeeman splitting will actually select (g2)tt as %n ground state. The 2n — > 2n + l 
ground-state-to-ground-state transition is then (g2)tt — Kti)t> as indicated left of the 
figure. It is the first line of a triple marking transitions to the lowest lying In + 1 states. 
Upon introducing a Zeeman energy, the spin-degeneracies of those are lifted, but only 
the two excited lines split in "V"-like manner, while the ground state-to-ground state 
transition (g2)tt~Kti)t has only one rightwards slanted branch (i.e. raises in energy 
with increasing field). The reason is that (g2)tt * s connected to the \c^) (= |t?4it)) ~ 
component of the energetically favored state (ti)-|- = ^(kt) + l^t)) + l c t) by in-tunneling 
of 4 - electrons. There is no possibility for a transition with f - electrons, hence a left 
branch does not exist. 

At a Zeeman splitting of y/2t/2, the process (g2)tt ~^ I t>t>t) becomes the ground 
state-to-ground state transition. The crossover is marked with (P) in Fig. [9^- Due to 
the in-tunneling of t - electrons, this resonance is continuously lowered in bias upon 
increasing the Zeeman energy further. 

At a Zeeman splitting of about 0.4 meV, we are exactly at resonance. As seen in the 
middle plot of Fig. [HJ a line triple has clearly separated from the ground state transition 
line. Upon comparison with Fig. [9^ it is immediately understood that it belongs to J, - 
electron transitions to the lowest lying 2n + 1 states. Concerning the corresponding |" 
- electron transitions, something interesting happens: In the point (P'), the left branch 
of the second "V" -shaped pattern, which belongs to the 2n + 1 state {tz)^, ends. The 
reason is that (t 3 )| consists of the same components as (ti)-(-. By in-tunneling of t - 
electrons, it can thus not be connected to (g2)th but rather to |g3) [= + \b))], see 

Fig. El sketch (P'). This 2n state is, compared to (g2)tt> lifted by the Zeeman energy 
and can only be populated by back-transitions from (ti)-(-. The state (ti)-j-, however, 
is not available below the transition (g2)tt ~~ (^i)t Fig- El sketch (P'), dashed arrow]. 
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(a) unpolarized (b) II polarized 

leads, P=0.0 leads, P=0.8 pjff. cond. 




Bias voltage eV bus [meV] Bias voltage ev^JmeV] 



Figure 9. Differential conductance vs. bias voltage and Zeeman splitting at a fixed 
gate voltage of approx. 2.0 meV, for an ACN transport setup with (a) unpolarized and 
(b) parallel polarized contacts. The sketch (P') on the left explains why the channel 
(g'2)tt ~ ^ (ti)t (dashed arrow) must be open in order to see transitions |g3)o (ts)-|-. 

This explains why the point (P') is positioned at the crossing with the resonance line 
marking this transition. 

Going on to a value of 0.8 meV of the Zeeman splitting, where the rightmost plot in Fig. 

E] is taken, we reside at low bias voltages within the 2n + 1 Coulomb blockade diamond; 

t 

the ground state transition is now the out-tunneling process |1\1\T) ~~ > (g2)th thus 
raising in bias for an increasing magnetic field. The behavior reverts again in the point 
(P"), where the Zeeman splitting has lowered the excited 2n state |t>tt>t) enough to 
change the ground state transition to 1 1, t> t)~ H t> TT> t): which involves out-tunneling 
of 4- - electrons. 

Finally, the Fig. [9b shows the data obtained if the calculation yielding Fig. [9^ is 
performed for ferromagnetic leads, polarized in parallel to the applied field. The only 
thing changing is the intensity of the lines. In particular, several of them are transformed 
into negative differential conductance lines. Such an effect is expected for any type 
of single electron transistor with parallel polarized contacts in magnetic field: upon 
opening a channel to a state from which the system can only escape via a weak (in our 
case J, - mediated) transition, NDC occurs as such slow processes hinder the current 
flow. The only exception are the ground state-to-ground state transitions, i.e. the 
edges of the Coulomb diamonds, where current starts to flow: to those, of course 
always positive differential conductance (PDC) lines belong. An obvious example is 
the |T, T, t) |1\T1\T) transition, which turns from NDC to PDC beyond the point 
(P"). 
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7. Conclusion 

We have studied the transport characteristics of fully interacting graphene armchair 
nanoribbons (ACNs) attached to ferromagnetic contacts. Short-ranged Coulomb in- 
teractions play an essential role in such systems, leading to an entanglement between 
bulk states and the ones localized at the zig-zag ends of the ribbons, thereby lifting de- 
generacies between various states. Importantly, the entanglement breaks the otherwise 
strict conservation of the bulk spin-S'z component, which leaves strong fingerprints in 
transport. 

The stability diagrams predicted for ACNs possess a two-electron periodicity and show 
already for a completely symmetric, unpolarized setup in zero magnetic field unique 
features like a characteristic transition line triple and pronounced negative differential 
conductance [29] . 

These effects, originating from the interaction-induced lifting and formation of states 
symmetric or anti-symmetric under exchange of the ribbon ends, are preserved for a 
parallel contact polarization. For an anti-parallel contact polarization, absence of vari- 
ous transition lines is observed due to spin-blockade effects and also the NDC features 
have vanished. The reason is that transition channels feeding the trapping state are 
disfavored, which leads even to a negative tunneling magneto-resistance. 
We have further investigated the transport behavior in magnetic field, for unpolarized 
as well as for in parallel polarized contacts. A change of the odd filling ground state with 
S z = h/2 to one with S z = 3H/2 is observed at a Zeeman splitting of \/2i/2, such that 
the value of the end-bulk exchange coupling t can directly be read off. Upon imposing 
a parallel contact magnetization, at several transition lines the differential conductance 
changes from the positive to the negative regime, because all 4- - mediated transport 
channels become weak. 

All in all, we found that short-ranged Coulomb interactions yield a strong influence 
of localized end states on the properties of ACNs. In particular, exchange makes the 
isolated bulk and end spin-S^ components a bad quantum number: only the sum of both 
is a conserved quantity. Due to this fact, ACNs might not be as ideal candidates for 
certain spintronic devices as previously regarded. On the other hand, the entanglement 
is a rich source of ACN specific features in transport. Recent achievements in fabrication 
of carbon nanostripes with defined geometries [151 [T9] raise the hope of an experimental 
confirmation of our predictions within the near future. 
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Appendix A. Diagonalization of H Q + V pp 



We start by rewriting Hq + V pp in terms of collective bosonic excitations. Concretely, 
we Fourier-expand the ID electron densities p ra (y) = i>ltr(y)i> r <r{y) > 

Mi/) = 2rE eiPW ^' ( A - x ) 

where the summation is over the wave numbers q = j^n q , n q £ Z. Then as shown e.g. 
in [12] the operators 

Kg ■= -^=Po q , & := ——p a -g, q>0, (A.2) 
^rT q q jrT q 

fulfil the canonical bosonic commutation relations [b aq , b\, i\ = S aa i8 qq /. 
As well known (12] the bosonization of Hq yields 



ft* = ^o (EE "Ata + \ E ^ ) > ( A - 3 ) 

V (7 q>0 cr / 



where A/" CT = Y~L c£„ c CTK „ counts the number of electrons with spin a. The first 
term in Eq. (1A.3I) accounts for collective particle-hole excitations, whereas the second 
term is due to Pauli's principle and describes the energy cost for filling up the spin 
degenerate single-electron states. Terms proportional to the total number of electrons 
have been omitted since they merely lead to a shift of the chemical potential in transport 
experiments. 

Bosonization of V pp ~ h can be achieved by rewriting the involved potentials in terms of 
electron densities and inserting the Fourier expansion Eq. flA.lj) . thereby making use of 
the definition Eq. (IA~2l) . We obt am 

T/b-b _ T/b-b , T/b-b , yh-b , T/b-b , T>b-b 
v oo — v f+ f+ + v f+ f- + v f- f+ + v f- f- + v hf+ 



rb—h t rh— b i T/b— b , T/b— b i t /b— b i T/b— b 

PP ~ V f+f+ + V/- + V f~f+ + V/" + V bf+ 



o-o-' g 

- I E E N (M* + + ^oA£ - ^ AC 2 , (A.4) 

with PVg as given in Eq. fl30|) . 

The last line of Eq. (|A.4[) describes the contribution of V pp ~ h to the system energy 
depending on the number of electrons in the two spin-bands. Here, E c = Wq is the ACN 
charging energy; J\ c = Jx^ + Jx± counts the total number of electrons. Spin alignment 
of electrons is favored by the term proportional to it, because it counteracts the energy 
cost for the shell filling in Eq. AA.3J) . 

Finally, the bosonized form of the end-bulk contribution Eq. (l2~5l) to V pp is 

V;; h = \[dy t p (y) E E v 7 ^ + h U) cos (^), ( A - 5 ) 

J a g>0 

which is linear in the bosonic operators, while those appear quadratically in Eqs. (1A.3j) 
and (EOj) . 
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In fact, any term of the form Eq. (IA.50 takes can be absorbed in the quadratic part 
of the Hamiltonian without any relevant impact on the spectrum, and we remain 
with H + V£~ h , which can be diagonalized in a standard way [12] by a Bogoliubov 
transformation [33]: One introduces new bosonic operators aj q and a^ q which relate to 
the old bosonic operators, b aq via 

b ag = K ( B j^ji + D ^) q ) q > 0, (A.6) 

K-± K-^. (A,) 

The transformation coefficients Bj q and Dj q can be expressed in terms of Xj q and Aj q , 
which were introduced in Sec. 13.31 

r> £ jq Xjq _ _ Aj q 



n r, ' 3i ~ I 

, £ i<? + Xjq) 2 — A q J (Ej q + Xj q ) 2 - A 2 

With our values, Eq. (1291) . for Xj q and Aj q we find approximately 



e cq = e n q ^Jl + W q /e , e sq = e n q \J 1 - (u/2e ) 2 « e n q , 
and for the transformation coefficients to the spin mode 

B sq = 1, D sq = 0. (A.8) 



The transformation coefficients for the charge modes depend, as for SWCNTs [351 HO] , 
on the ratio g q = eo q /e cq : 

B c q = \{V^+ Vv^) > A, = \{y/9~ q - • ( A -9) 

Exploiting these relations yields then the diagonalized Hamiltonian Eq. (J2BD- 

Appendix B. The matrix elements of the non-diagonal bulk-bulk interaction 

The evaluation of the non-diagonal bulk-bulk terms follows closely the procedure applied 
for SWCNTs [33]. In order to calculate the matrix elements of the non-diagonal bulk- 
bulk interaction we have to derive an expression for the matrix element 



Na e m 



' npp 



N'a e, m') =: M [r]W/ _ (N, a e , m, N', a e ' , m', y) 



Na e m! 



■= 5jv j7 v/ S#e,#«> M [r ] [a ] f _(N,y)M[ r ][ a ] f _(m,m',y), 

which we have factorized in the last step into a fermionic and a bosonic part. We 
express the operators ipraiy) in terms of the bosonic operators b aq and b\ , q > 0, using 
the bosonization identity |42j . 

AMv) = f}aK ra {y)^ Uv) ^ ra{v) - (B.l) 
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The operator fj a is the so called Klein factor, which annihilates an electron in the 
cr-branch and thereby takes care of the right sign as required from the fermionic anti- 
commutation relations; in detail, 



N - e a ,m) . 



K rrT (y) yields a phase factor depending on the number of electrons of spin a, 



K Ta {y) 



1 i£ S gn(r)(.X?- .+i)j/ 



Finally, we have the boson fields i<p r a{y), 

The fermionic part is then given by 

and the bosonic part reads 

Mm(fh,fh',y) = (m| e~ 4 h M e -tin<v) e -$i a <v) x 



M {l] (N, y ) = (n iqM€K(v)&K*(v)fihKu(y)fh 



N' 



(B.2) 



(B.3) 



(B.4) 



(B.5) 



(B.6) 



In order to improve readability we have replaced the indices ro~ by a single index I. 
Using the relation flB.2j) for the Klein factors fj a , together with the fact that S a — f~, 
[o~]f- = [a, — a, — a, a], and the definition Eq. ( 1B.3|) of the phase factor K ra {y), it is 
straightforward to show that 

1 



Mr, 



] f -(N,y) 



(2Ly) 



>QN[r]a(y)i 



(B.7) 



where Qtf [r]a (y) = exp 
Hence, for S r = u, [r] u = 



y 

r, r 



if^ yN a sgn(r 4 - r x ) - AL CT sgn(r 3 
r, — r] , we obtain 



„ \ i sgn(r4+r 3 -r2-ri) 

r2) H 2 



Qn[t)uM = eX P 



i^sgn(r) (N <r + N- (r + l)y 



it 



which is oscillating fast with N c = + N± and thus completely suppresses the S r 
contribution away from half-filling. The only remaining term in V^~ p h is consequently 
S r S a = bf~, for which we get with [r] b = [r 



QN[r] h a(y) = eX P 



2tt 



-i— sgn(r) (N a - N y 



We can now restrict our further analysis to the bosonic part (fn, ml, y). We 

are going to express the fields i4> r a{y) i n Eq- (IB.6j) in terms of the bosonic operators 



a 



and subsequent normal ordering, i.e., commuting all annihilation operators a 
to the right side and all creation operators a},- q to the left side. In a first step we use 
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the relation e ^ (j/) e^ t(i/) = e^^e^'We^^^, following from the Baker-Hausdorff 
formula, to obtain from Eq. (IB.6|) 

/ ~ 4 ~l ~ 4 

g-*En=:l>I„ (!/) g-iE„=l<fon (V) 



M[q(m,m',y) = C w (y) (m 



7? I 



(B.8) 



^,4 



where ^ =1 </>z„ denotes the sum 0^ + 



and 



q (y) — e I*^3 (f )'^« T 4 (f )1 e [ _ ^2 (y)'^? 3 fe)+^? 4 (f)] gh^ii (y)- i( t>J 2 {y)+ i 4>l i iv)+ i 4>l 4 (y)] 



With the definition Eq. (IB .41) of the boson fields, we can easily derive the anti- 
commutator relation 



i4>ra{y),i4>t> a >(y') = -8a,a>^2 



Jq(ry-r'y') 



/n, 



(B.9) 



<3>0 



which allows us to simplify C [r]b[a]f _ (y) = exp [-i}_ ra (y), i<j>l a (y)] exp [-i(p ra (y), i$_ ra {y)] 
(1 - exp [-2irTxylL y \Y l (1 - exp [2irny / L y \y l = [4 sin 2 {wy/L v )]- x . 



Applying the Baker-Hausdorff formula once more, we obtain further for the second 
contribution to Eq. (IB. 8 1) 



-*En=i# n (v) e -*E„ =1 <k (y) = e -iE„=i (A„ &/)+#,, (v)) J (v),iE»'=i& n , (») 



Using the definition Eq. (1B.4h together with the transformation between the operators 
Oq-(j and d jg , Eq. (1A.6[) . we get 



l jqj ■ 



(B.10) 



jq>0 



In terms of A J a , Bj q and Dj q , which were introduced in Appendix A the coefficients 
Xii(x) read 



\H(y) = ( e isen(r)qy B jq - e- isgn( - r)qy D jq ) 

\ /Tin 



(b.ii; 



and plugging in the corresponding values, cf. Eqs. (1A.7[) . ( 1A.8j) and (IA.9j) . it is easy to 
calculate 

~ x >) ■= -EL W (B.i2) 

for [I] = [r] b [a] f -. We find that 



x ikw] f - (y) = ' a m 6 h / - (f) = ~ 2i \ — s s n ( r(T ) t^+&J sin (^)- 



1 ssO 

-4 



Again using the Baker-Hausdorff formula yields for the exponential e i En=i(^n(s/)+^„(s/)) 
= exp (-1 £ 3 >ol4 9 ]6W/ _ ( y )|») e -^>°%v ^e^-^w,- W*- > such that in total 



~ 4 -| - 4 

e -*E„ = l</'r„ CT „(?/) e -*E„ = l</ > r nC T n (»/) 



771 



5 



(B.13) 
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An explicit evaluation shows 



r , (y) 



A [r]b[t7]f _(y) = e E -^( 2 - e " 2 ^- e2!r9!/ )e- 4E -^ sm2(w) = 1. 



The function F(X,m sq ,mf s 



sqj 



(rh s 



F(X, m, m!) = \ Q(mf - m)X m '- m + Q(m - m! 

(-IAI 2 )' 



m' s ) is given by |1QI 135] 



X 



m 



Or 



(B.14) 



where m m i n / max = min / max(m,m'). Combining Eq. ( IB. 81) and Eq. ( IB. 131) we finally 
obtain 

Sri 



Mr 



J m c ,m' c 



4 sin (ny/^ 



(B.15) 



Altogether, we get with Eqs. (1B.7|) and (1B.15|) to an expression for the matrix elements 

rb-b 

npp 



of Vn PP h away from half-filling 



(Na e m 



b-b 
npp 



a r/ -»e/ at / 

Nam 



(Na e rh 



V bf- 



N'a e 'm' 



ii 



2L„ 



■s^s^s^, [ d y j^jz^z ) II F C x M b[ * ]f - (y)> m ^ < 



For ^ q \m sq — m' sq \ < 1, the evaluation of this expression is problematic as 
the divergence arising from l/[4 sin 2 (ny/Ly)] remains uncompensated. Hence, the 
evaluation of the corresponding matrix elements needs special care. The origin of this 
divergence lies in the fact that, if no bosonic excitations are present, the N conserving 
processes depend on the total number of electrons in the single branches [compare to 
the fermionic contributions to Hq + V pp in Eq. (128j) ]. Since the bosonization approach 
requires the assumption of an infinitely deep Fermi sea [12], this leads, without the 
correct regularization, necessarily to divergences. These findings are in complete analogy 
to the theory for SWCNTs [35]. In the following we exemplify the proper calculation 



for ( Nrh 



V, 



b-b 



Nmj. 



Regularization of M^.] b [ a ] _ (N, a e , rh, N', a e> , rh', y) for rh — rh' Regularization for the 
matrix elements of the non-density-density bulk-bulk interaction is needed in case of 
J2j q \ m jq — m ' jq\ < 2, since in that situation M^ b [ a ] f _(N,a e ,rh, N',a e ' ,fh' ,y) diverges 
due to the factor 1/4 sin 2 (j^-y) in Eq. (1B.15j) . 

Here we show the details of the proper regularization for rh = rh'. In this case we make 
the expansion 

J] F(AJ5 bW _(y), rh, rh) = 1 + G(sin 2 ), 
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where C(sin 2 ) contains only terms (sin(gy))' 9 with ^2 q t q > 2 and which 'cure' 
the 1 / sm 2 (%y / L y ) divergence appearing in M^ b ^ _(N,m,m,y). Therefore we are, 
compare to Eq. ( IB. 161) . left with the regularization of 



g Ly g Ly 



4 sin 2 (iry/L y ) 4 sin 2 (ny/L y ) 1 _ jr%v 1 _ e -*%v 



with the second equality obtained following [35] . Using further ^rL-oo e " 
this leads to 



-vnx e 



QN[r] b a(y) 

4 sin 2 {ny/Ly 



E 



in'r^-y 
V > g y . 



E 



Integration over y brings us to 

QN[r] b a(y) 



dy 



4 sin 2 (ny/Ly 



n=— oo n'=— oo 



Summarizing, the regularized expression reads: 



( Na e m \v£- h \ N'a e m )=u5^ 



J]min(A^ (T ,iV_ (T ) 



+ 4E/^j§^n(^ M/ _ M , m .,, m , 9 )-i 



Appendix C. The matrix element of the non-diagonal end-bulk interaction 



Also for the non-density-density end-bulk-scattering, we omitted the calculations in 
the main part of the text and give the detailed evaluations here. We have to 
evaluate matrix elements of Eq. ( 1261) . of the form M p arla ,(N, a e , rh, N', a el , rh') :- 
(^Na e fh 
previous 
end part 



^iAyp)^r>v>(y P )Y.K X di lpKx dv P K X N'a e 'fh'j, with y_=0, y+=L y as employed 
y. We can factorize the matrix elements into a fermionic, a bosonic, and an 



ilC rV (TV, a e , rh, N>, a e ', rh') = (TV, N')Ml r , a , (rh, rh')M p aa! (a e , 3") 

The easiest to give is the end contribution: 



a 



dl 



a 



6 a <L a <LI Sfje a l 6 a el 



(C.l) 



The end electron operators act on the configuration at the p-end only, trying to transform 
the spin from a to a', while the spin at p — —p must be untouched. 

What remains is (iVm|^J CT (|/p)^ r / (T /(?/p)|iV / m'). Upon applying the bosonization 
identity Eq. (IB. II) to rewrite ip rcr (y p ), it becomes a product of the fermionic part 
M p ar , a ,(N, N'), which is calculated below, and the bosonic part, 



M p ,(rh,rh') = (rh | e -^fe) e -^ ife) e^ T ' fe) e^' fe) 

= C w (y p )A w (y) JJ F(\%(y),m jq ,m' jq ), 
j,g>o 



111 
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where we abbreviated ra by I and r'o' by V and exploited Eq. (IB. 9ft to obtain Cw(y) : = 

e [-4i(y),i4>l,(y)} _ e *<r,a' E 9>0 ^ i(r " r ' ) " = e ^,a'E g >o^-{^(25 r> _ r , S 2/)+i S gn(r)sin(2a r , > _ r ., g2 /)}_ 



the remaining exponentials we did as for the bulk-bulk scattering in |Appendix B 



and while F(\,m,m') is the function defined in Appendix B Eq. ( 1B.14I) . Au>(y) : = 
e IK^W-^W),<^W-^W)] e -|E 3 ,|^| 2 = e E 9 >o^(i-^ iCT ,cos(25 rir , gy )) e _i Ej jA^( 3/ )| 2 guch 

that 

A u ,{y)C lv {y) = e E ^ i^ 1+ ^-' sgn(r)sin(2w) )e-^-l^r. (C.2) 

Further, we set in Eq. (IB. lip to expand A^(y) := —\j 9 (y) + \\?{y) = 
-|= [B jq (AJ, e ^gn(r')w _ e isgn(r)w ) - D jq (A^ e - issn(r ' )(ra - A^e- isgn(r)w )] . For our 
matrix element M[ l ,(m,m') we have only to consider the two cases y = y_ = and 
y = y + = L y , where the expression simplifies under application of Eqs. (1A.7[) and (lA.8h 
along with sin(O) = sin(±gLj,) = 0, cos(0) = 1, cos(±qL y ) = —1 to 



Hi'iVp) = $j,s8<T,-a>sgn(pcr)\/2/n q . 
Using this result in Eq. (IC.2p yields 



A w {y p )C w (y p 



const. a' = a, 
1 a' = -a. 



With this ingredients we arrive eventually at the crucial result 

, ,„ , _ —* r\ const. a' = a, . . 

MfjAm.m ) = < . „ _,. . . , — , — , , , (^•3) 

| $rh c ,m' c U q>0 ^(sgn(pa) y/2/n q , m sq , m' sq ) a 1 = -a. 

Finally, the fermionic part is given by 



N 1 

where the Klein factors fj a , Eq. (1B.2j) . and the phase factors K ra (y), Eq. (1B.3|) . were 
straightforward to evaluate. We split now the sum contained in V°~p , Eq. ( 1261) . 

and therewith the further analysis, into two cases. 
First, let a' = a hold. Then, 



TUT U \ ' / 



6 n,n'/( 2L v) y P = y~ = o, 

~ rr ' S N,N'/ ( 2L v) Vp = y+ = L v 
The total contribution related to this case is (t = i np ) 

According to Eq. ( 1C3[) . in fact M^ r , a (m, m!) = M~ ar , a (m,rh') = const, and so this 
contribution identically vanishes. 
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We are left with the part of the sum where a' a, for which we can read off the following 
fermionic contributions from Eq. flC.4j) : 

6, 



M? ar ,_ a (N,N>) 



1L,, 



-sgn(a) - 6 p , + rr'(-lfl) , (C.5) 



To obtain the first equality we used (—1)^ <5<t4 = — (— l) Ar i'sgn(o"), while for the second 
one we needed exp {m[r'(N'_ a + 0.5) - r(N' a + 1.5)]} = -sgn(rr')(-l) JV '»(-l) JV '-'. 
Putting Eqs. ( IC.lj) . ( 1C.3j) and ( 1C.5I) together, the final result for the non-diagonal 
end-bulk-scattering is: 



N a e = rcr_i_ , <r_) rh 



npp 



N' a e ' = (a'a'_) m' 



-l) N r - cL 



■l) N i 



sgn(a p )l[F (psgn(a' p )J2/ 



"n q ,m sq ,m sq 
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